The Adiabatic Instability on Cosmology's Dark Side 
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We consider theories with a nontrivial couphng between the matter and dark energy sectors. We 
describe a smah scale instability that can occur in such models when the coupling is strong compared 
to gravity, generalizing and correcting earlier treatments. The instability is characterized by a 
negative sound speed squared of an effective coupled dark matter/dark energy fluid. Our results are 
general, and applicable to a wide class of coupled models and provide a powerful, redshift-dependent 
tool, complementary to other constraints, with which to rule many of them out. A detailed analysis 
and applications to a range of models are presented in a longer companion paper. 

PACS numbers: 



In order for our cosmological models to provide an ac- 
curate fit to current observational data, it is necessary to 
postulate two dramatic augmentations of the assumption 
of baryonic matter interacting gravitationally through 
Einstein's equations - dark matter and dark energy. A 
logical possibility is that these dark sectors interact with 
each other or with the normal matter [H, 0, Q • ^ number 
of models have been proposed that exploit this possibility 
to address, for example, the coincidence problem. 

Such models face a range of existing constraints arising 
from both particle physics and gravity. In this letter we 
consider perturbations around the cosmological solution 
and demonstrate the existence of a dynamical instability 
which we term the adiabatic instability. This instability 
is characterized by a negative sound speed squared of the 
effective coupled fluid P, \B\ and was first discovered Q 
in a context slightly different to that considered here - 
the mass varying neutrino model of dark energy. Our 
aim here is to give a general treatment of the instabil- 
ity, applicable to a wide class of models, to identify the 
regimes in which the instability occurs, and to delineate 
the resulting redshift-dependent constraints. 

Class of Models: We begin from the following action 
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where g^v is the Einstein frame metric, </> is a scalar field 
which acts as dark energy, and are the matter fields. 
Here we have adopted a signature (-,+,-f ) and defined 
the reduced Planck mass by ra^ = (SttG)^^. The func- 
tions aj(0) are couplings to the matter sector. This 
general action encapsulates many models studied in the 
literature 0. The field equations are: 
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[(Pj +Pj)Uit.u-iu+Pi9i,u] ,(2) 



where we have treated the matter field(s) in the j sector 
as a fluid with density pj and pressure pj as measured in 
the frame e^"-i g^y , and with 4- velocity uj ^ normalized 
according to g^'^Uj^Uji, = — 1. 

We consider models with a baryonic sector (ah(0)) and 
a composite dark matter sector, with one coupled species 
with density and coupling adfj)) = ^ind another 

uncoupled species with density pco and coupling aco — 
0. Neglect the gravitational effect of the baryons, using 
pc = pco = 0, and defining pj = e^"' pj gives 



(4) 



and VpV^ 



= 0, where we have defined an 
effective potential by V;ff((/)) = + e"('^)/9c • The fluid 
obeys V^(pcO = 0, and <V^< = -(5^'' + u>c)Vi.a. 

Adiabatic regime: The effective potential Kjff((/)) may 
have a minimum resulting from the competition between 
the two distinct terms. If the timescale or lengthscale 
for (f) to adjust to the changing position of the minimum 
of Voff is shorter than that over which the background 
density changes, the field (f> will adiabatically track this 
minimum Q. In this case the coupled CDM component 
together with (p together act as a single fluid with an 
effective energy density pcft and effective pressure PcS'- 



PcsiPc) 



ea[0„,(Po)]^, 
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PesiPc) = -V[(f>,n{pc)] . (6) 

Here 4>ai{pc) is the solution of the algebraic equation 

K'ff(0) = X/'(0) + a'((/))e"(^)pe = O (7) 

for (f). Eliminating pc between Eqs. ^ and ^ gives the 
equation of state poff = Poff (poff)- 
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For cosmological background solutions, we assume that 
the coupled fluid acts as the source of the cosmic accel- 
eration. In the adiabatic approximation, the effective 
fluid description is valid for the background cosmology 
and for linear and nonlinear perturbations. Therefore, 
the equation of state of perturbations is the same as that 
of the background cosmology, and the matter and scalar 
field evolve as one effective fluid, obeying the usual fluid 
equations of motion with the given effective equation of 
state. 

A necessary condition for the validity of the adiabatic 
approximation is that the lengthscales or timescales £ 
over which the density pc varies are large compared to 
inverse of the effective mass 
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of the scalar field. More precisely, we can show that the 
condition is [sl 
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this condition is necessary to justify dropping the terms 
involving the gradient of cf) from the fluid and Einstein 
equations. In most situations the logarithmic derivative 
factor is of order unity and can be neglected. In Rcf. Q 
we also derive a non-local sufficient condition for the va- 
lidity of the approximation, which generalizes conditions 
in the literature for the chameleon (thin-shell condition) 
and f{R) modified gravity [l2l | models. Condition 
(PI is not very stringent; many dark energy models admit 
regimes where it is satisfied for the background and for 
linearized perturbations over a range of scales. 

In the adiabatic regime, the inferred dark energy equa- 
tion of state parameter in the case as = is 
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with ofQ = ct{4>Q) the value today. Thus, w is precisel 
today, and gcnerically satisfies u; < — 1 in the past 0, Q ■ 

Adiabatic instability: We write the potential V{(f) as a 
function V{a) of the coupling function a{(j)) by eliminat- 
ing 0. This gives, from Eqs. ^ and (O, 
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dV/d(t) _^ dV 
da I d(\) da 
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The square of the adiabatic sound speed, = P/p is 
then given by 
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dpeff _ dpcs/da 
dpcs dpcs/da 
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dot 
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In the adiabatic regime the effective sound speed, re- 
lating to local perturbations in pressure and density. 



c^(fc, a) = 6P{k,a)/Sp{k,a), tends towards the adiabatic 
sound speed and is always negative, since dV/ da must be 
negative so that Eq. ([7]) admits a solution, and d'^V/da'^ 
must be positive so that ([5]) yields a positive m^g. From 
here in, we consider the regime in which this adiabatic 
limit has been reached, and take = c^. 

Consider now a perturbation with lengthscale £. In 
order to be in the adiabatic regime we require £ ^ . 
The negative sound speed squared will cause an exponen- 
tial growth of the mode, as long as the growth timescale 
~ £/ \/|Cs| is short compared to the local gravitational 
timescale nip/ y/pcsipc)- Combining Eqs. (O, ([7]) and 
® yields c>2^ = (a')'K« = ia')^Pes/{V/Vc. - 1), 
and therefore the instability will operate in the range 
of lengthscales given by 



rricffipc) 
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Here the quantity dhiV/da{a) on the right hand side is 
expressed as a function of cj) using a = a{(f)), and then as 
a function of the density using = (j),^{pc). In order for 
this range of scales to be non empty, the dimensionless 
couphng mpla'l must be large compared to unity, i.e., 
the scalar mediated interaction between the dark matter 
particles must be strong compared to gravity. 

There are two different ways of describing and under- 
standing the instability, depending on whether one thinks 
of the scalar-field mediated forces as "gravitational" or 
"pressure" forces. In the Einstein frame, the instability is 
independent of gravity, since it is present even when the 
metric perturbation due to the fluid can be neglected. In 
the adiabatic regime the acceleration due the scalar field 
is a gradient of a local function of the density, which can 
be thought of as a pressure. The net effect of the scalar 
interaction is to give a contribution to the specific en- 
thalpy h{pc) = J dp/ Pc of any fluid which is independent 
of the composition of the fluid. If the net sound speed 
squared of the fluid is negative, then there exists an insta- 
bility in accord with our usual hydrodynamic intuition. 

In the Jordan frame description, however, the insta- 
bility involves gravity. The effective Newton's constant 
describing the interaction of dark matter with itself is 
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where k is a spatial wavevector [8|. At long lengthscales 
the scalar interaction is suppressed and Gcc ~ G. At 
short lengthscales, the scalar field is effectively mass- 
less and Gcc asymptotes to a constant. However, when 
77ip|a'| S> 1 there is an intermediate range 
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over which the effective Newton's constant increases like 
Gcc oc k^. This interaction behaves just like a (negative) 
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pressure in the hydrodynamic equations. This explains 
why the the effect of the scalar interaction can be thought 
of as either pressure or gravity in the range of scales 
Note that the range of scales ([T5|) coincides with with the 
range (jl3p derived above, up to a logarithmic correction 
factor. 

From this second, Jordan-frame point of view, the in- 
stability is simply a Jeans instability. In a cosmological 
background the CDM fractional density perturbation tra- 
ditionally exhibits power-law growth on subhorizon scales 
because Hubble damping competes with the exponen- 
tial (Jeans) instability one might expect on a timescale 
of Xj \JGp. In our case, however, the gravitational self- 
interaction of the mode is governed by Gcc{k) instead of 
G, and consequently in the range (fT5|) where Gcc ^ G 
the timescale for the Jeans instability is much shorter 
than the Hubble damping time. Therefore the Hubble 
damping is ineffective and the Jeans instability causes 
approximate exponential growth. 

Examples of Models: For single component dark matter 
models, one can find coupled models in the adiabatic 
regime [3, fl^ . However in the strong coupling limit 
mp|a'| 3> 1 of interest here, they typically do not yield 
acceptable background cosmologies. Therefore we focus 
on composite dark matter models. 

As a first example we consider a constant coupling 
function and an exponential potential 

a(,/)) = -/3C-^ , V = Voe-^'''/'"^ , (16) 

TOp 

where /3 = •\/2/3 and C < and A are constants. 
The Friedmann equation in the adiabatic limit is then 
3mpH^ = V + e"pco + Pc, in which the first two 
terms on the right hand side act like a fluid that, for 
\C\ ^ 1, approaches a cosmological constant. Thus, 
the background cosmology is close to ACDM for large 
enough |C|. Since the fraction of coupled dark matter 
is flco = e°'pco/{3mpH'^), in the asymptotic adiabatic 
regime, ilco = A(l — ilc)/{\ — f3C), and flc ^ 0.3 to- 
day, ilco must be small for large coupling, |C| 1. If 
the parameters of the model are chosen so that fie ~ 1 
today, then the maximum and minimum lengthscales for 
the instability are £max ~ and >Cmin ~ {HoP\C\)~^ . 
Taking the Jeans view, it is then possible to show [8[ that 
the instability should operate whenever modes arc inside 
the horizon and in this range. 

These expectations are confirmed (figure [1} by a nu- 
merical analysis of a two component coupled model. We 
use A = 2, C = -20, Hq = 70kms-iMpc~\ baryon 
fractional energy density, fib = 0.05, uncoupled CDM 
component, flc = 0.2, coupled component, flco = 0.05, 
and potential fractional energy density, fiy = 0.7. We 
fix initial conditions of (p/m-p = 10^^" and = at a = 
10""'^° (initial conditions at least within (j>/mp = 10~^° — 1 
give the same evolution because of a scalar dynamical at- 
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FIG. 1: [Bottom] The two component coupled dark energy 
(CDE) model, with A = 2 and coupling C — —20 with 
Ho = 70kms-iMpc-\ = 0.05, = 0.2, Q^o = 0.05, 
and Slv = 0.70. At late times the scalar field finds the adia- 
batic minimum with asymptotic equation of state, and sound 
speed = —1/(1-1-7) = —0.89, able to reproduce a viable back- 
ground evolution consistent with supernovae, CMB angular 
diameter distance and BEN expansion history constraints. 
The figure shows the evolution of the effective equation of 
state, Weff = Ptot/ptot ~ {2/3){d\nt/dlna) — 1, (black full 
line), adiabatic speed of sound, — Ptot/ptot, (blue long 
dashed line) and effective speed of sound for = SPtot/Sptot 
at A: = 0.01/Mpc (red dot-dashed line). The effective equa- 
tion of state for a comparable ACDM model with fie = 0.25, 
Q,b ~ 0.05 and Qa ~ 0.7 is also shown (black dashed line). 
[Top] The growth of the fractional over-density S = Sp/p for 
k = 0.01/Afpc for the coupled CDM component, S^o, (red 
long dashed line) and uncoupled component, 5c, (black full 
line) in comparison to the growth for the ACDM model (black 
dashed line). At late times the adiabatic behavior triggers a 
dramatic increase in the rate of growth of both uncoupled and 
coupled components, leading to structure predictions incon- 
sistent with observations. 



tractor) and assume that the CDM components have the 
same initial fractional density perturbations Sc = Sco, 
fixed by the usual adiabatic initial conditions. As shown 
in the bottom panel of figure [1] the background evo- 
lution is consistent with a ACDM like scenario, with 
Wctf ~ —0.69 today, approaching Wcs ~ —0.89 asymp- 
totically. In the top panel we see that, once the scalar 
field has entered the adiabatic regime, giving rise to ac- 
celerative expansion, the density perturbations undergo 
significantly increased growth, in stark contrast to the 
ACDM scenario in which accelerative expansion is typi- 
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cally associated with late-time suppression of growth. 

In summary, these models provide a class of theories 
for which the background cosmology is compatible with 
observations, but which are ruled out by the adiabatic 
instability of the perturbations. 

Another interesting class is the chameleon models 0, 
[lo| for which the adiabatic regime has been previously 
demonstrated in static solutions for macroscopic bod- 
ies like the Earth, and also in cosmological models 11|. 



One well-studied example of these has inverse power law 
potentials, together with the constant coupling function 
in p6|) . for which the effective potential is then 
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,Pc) = XM' 



M 



(17) 



where M is a mass scale and ti > and A are constants. 
The existence of a local minimum, and hence an adiabatic 
regime, in (|17p requires C < 0. We shall restrict attention 
to the regime > pcrit = n\M^{^-(iCM jm^Y . The 
sound speed squared is 
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which is always negative as expected. 

The range of spatial scales £ over which the instability 
operates for a given density pc ^ Pcrit is non-empty for 
/3|C| » 1, and is given by 



1 < 



(n-f l)(/?C)2perit / P 



Pcrit 
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If 4> behaves as dark energy, we require Pcrit ^ -f^o™p- 
Then for pc ^ Pcnt , the maximum lengthscale is of order 
Hq^, and the minimum is ~ {HoP\C\)~^ . Thus a large 
set of cosmological models are in the unstable regime at 
Pc ~ /Ocrit (if P\C\ 3> 1), ruling them out in this regime. 

In this letter we have demonstrated the existence and 
broad applicability of the adiabatic instability - operat- 
ing in models in which there exists a nontrivial coupling 
between dark matter and dark energy. We have pre- 
sented general expressions for the conditions under which 
the adiabatic instability is relevant, and, when so, the 
lengthscales over which it operates. This work provides 
a new way to constrain interactions in the dark sector, 
and heavily restricts the class of models consistent with 
cosmic acceleration. 

In a companion paper Q , we derive in detail the results 
presented in this letter, and apply the results to a wide 
class of coupled models including couplings to both CDM 
and neutrinos. 
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